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Abstract 

The Lie group method provides an efficient tool to solve nonlinear partial differential equations. This paper 
suggests a fractional partner for fractional partial differential equations. A space-time fractional diffusion 

r^ , equation is used as an example to illustrate the effectiveness of the Lie group method. 
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I/"") ■ 1 Introduction 

>' 

In the last three decades, researchers have found fractional differential equations (FDEs) useful in 
various fields: rheology, quantitative biology, electrochemistry, scattering theory, diffusion, transport 
theory, probability potential theory and elasticity [1], for details, see the monographs of Kilbas et 
al. [2], Kiryakova [3], Lakshmikantham and Vatsala [4], Miller and Ross [5], and Podlubny [6]. On 
the other hand, finding accurate and efficient methods for solving FDEs has been an active research 
r ^\ ' undertaking. 

Since Sophus Lie's work on group analysis, more than 100 years ago, Lie group theory has become 
S^ more and more pervasive in its influence on other mathematical disciplines [7, 8]. Then a question 
may naturally arise: is there a fractional Lie group method for fractional differential equations? 

Up to now, only a few works can be found in the literature. For example, Buckwarand and Luchko 
derived scaling transformations [9] for the fractional diffusion equation in Riemann-Liouville sense 

d a u(x,t) ^d 2 u(x,t) ^ . . 

~^ 1 = D ^x^' 0<a,0<x,0<t,0<D. (1) 

Gazizov et al. find symmetry properties of fractional diffusion equations of Caputo derivative [10] 

d a u(x,t) , d(k(u)u x (x,t)) , . , 

g ta =k K V y > 0<a,0<x,0<t,0<k. (2) 

Djordjevic and Atanackovic [11] obtained some similarity solutions for the time-fractional heat 

diffusion 

d a T(x,t) , d 2 (T(x,t)) 

dta = k dx 2 — . 0<a,0<x,0<t. (3) 

In this study, we investigate anonymous diffusion [12] 

d a u(x,t) d 2l3 u(x,t) 



dt a dx 2 ? 



0<a,P<l,0<x,0<t, (4) 
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with a fractional Lie group method, and derive its classification of solutions. Here the fractional 
derivative is in the modified Riemann-Liouville sense [13] and — »"ffi ' is defined hy4^{ — ,wj' ' ) . 

2 Characteristic Method for Fractional Differential Equations 

Through this paper, we adopt the fractional derivative in modified Riemann-Liouville sense [13]. 
Firstly, we introduce some properties of the fractional calculus that we will use in this study. 

(I) Integration with respect to (dx) a (Lemma 2.1 of [14]) 

oOW = ^r [ X (x-O a - 1 f(m = V A- [ S I* f(®m a ,0<a<l. (5) 

T{a) J T(a + l)J 

(II) Some other useful formulas 

d f„Mu^ ,n«J_ l '('+^) T /3-a_ ^ 



dx r(l + p — a) 



The properties of Jumarie's derivative were summarized in [13]. The extension of Jumaire's 
fractional derivative and integral to variational approach of several variables is done by Almeida et 
al. [15]. Fractional variational interactional method is proposed for fractional differential equations 

[16]- 

It is well known that the method of characteristics has played a very important role in mathemat- 
ical physics. Preciously, the method of characteristics is used to solve the initial value problem for 
general first order. With the modified Riemann-Liouville derivative, Jumaire ever gave a Lagrange 
characteristic method [17]. We present a more generalized fractional method of characteristics and 
use it to solve linear fractional partial equations. 

Consider the following first order equation, 

a(x,t)^ + K X ,t)^ = c(x,t). (7) 

The goal of the method of characteristics is to change coordinates from (x, t) to a new coordinate 
system (xq, s) in which the PDE becomes an ordinary differential equation along certain curves in 
the x—t plane. The curves are called the characteristic curves. More generally, we consider to extend 
this method to linear space-time fractional differential equations 

a(x,t)^^ + b(x,t)^^- = c(x,t),0<a,/3<l. (8) 

With the fractional Taylor's series in two variables [13] 

Similarly, we derive the generalized characteristic curves 

(10) 

(11) 

(12) 
Eqs. (10)-(12) can be reduced to Jumarie's result if a = f3 in [17]. 



du 

d^ = c{x 


■;t), 




{dxf 


-a(x 


,t), 


r(i+p)ds 


(dt) a 


- b(x 


,t). 
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As an example, we consider the fractional equation 



' ""«(*■"+ /"V 8 "?*-"- 0,0 <a,f><l. (13) 



ft: 



T(l + /3) dxl 3 r(l + a) dt 
We can have the fractional scaling transformation 

x P 2t a 

" = u( F(I+^) / r(i+^ ) - (14) 

2 

Note that when a = f3 = 1, as is well known, ^- is one invariant of the line differential equation 

ox at 

3 Lie Group method for Fractional diffusion equation 

With the proposed fractional method of characteristics, now we can consider a fractional Lie 
Group method for the fractional diffusion equation, which are the generalizations of the classical 
diffusion equations treating the super-diffusive flow processes. These equations arise in continuous- 
time random walks, modeling of anomalous diffusive and sub-diffusive systems, unification of diffusion 
and wave propagation phenomenon [18 - 23]. 

We assume the one-parameter Lie group of transformations in (x, t, u) given by 

r(i+/?) = T(TTW) ^^ x ' + ^( e /i 

r(fe) = r(fe)+ er ( x ' t ^) + ( e )' (16) 

u = u + ecj)(x, t, u) + 0(e) , 

where e is the group parameter. 

We start from the set of fractional vector fields instead of using the one of integer order [9 - 11] 

V = £,(x,t,u)DP + T(x,t,u)D? + <t>(x,t,u)D u . (17) 

The fractional second order prolongation PA 2 "'Vofthe infinitesimal generators can be represented 
as 

dDfu 9 QD^u dD^u V QDfu dD^Dfu ^ ' 

As a result, we can have 

Pr^V(A[u}) = 1 (19) 

, A r i d a u(x,t) d^uixj) 

where A\u] = ^ > q^t 1 - 



Eq. (19) can be rewritten in the form 



(*M 



= 0. (20) 

Aful=0 



The generalized prolongation vector fields are defined as 

<t>® = D?<t>-(D?OD?u-(D?T)D?u, (21) 

<f>^ = D^-{D x e)D^u-{D^r)Dtu, (22) 

^H = D^-2(Df^)D^-(D^)Dfn-2(Dfr)DfD>-(Dfr)D> t . (23) 

Substituting Eqs. (21)-(23) into Eq. (20) and setting the coefficients to zero, we can obtain some 
line fractional equations from which we can derive 
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€(x,t,u)-a + c 4 r( i +/3) + 2c 5 r(1+a) + 4c 6 r(1+ ^ r(1+a ) , 
r(a,t,M) = C2 + 2c 4 r^ 5 y + 4c 6 r( ff 2a) , 

0(x,t,n) = (c 3 -c 5 r(f T? y-2c 6r ^y-c 6 r(1 ^ 2)g) )n + a(x,t), 

where Cj(z = 0...6) are real constants and the function a(x,t) satisfies 

d a a(x,t) d 2/3 a(x,t) 



dt a dx 2 ? 



, 0<a<l,0</3<l. (24) 



It is easy to check that the two vector fields {Vi, V2, Vj$, V4, V5, V s } are closed under the Lie bracket. 
Thus, a basis for the Lie algebra is {Vi, V^Vj, V4, V5}, which consists of the four-dimensional sub- 
algebra Oi,V2,V3,Vk} 



„ — d? „, _ a a „, _ a „, _ x^ a' 3 I 2t a 9" 
Ul ~~ a^ 3 " ' y2 — at a ' U3 ~~ an ' U4 — r(i+/3) &? + r(i+a) at« ' 

_ 2t a a' 3 mC g 

u 5 — r(i+«) Th? ~ r(i+/3) a^' 

_ 4t a x* 3 a' 3 4< 2c * a™ / x 2/3 , 2t a \ a 
U6 — r(i+«) r(i+/3) a^ + r(i+2«) at^ ~~ vr(i+2/3) + r(i+«) J n a«"' 



and one infinite-dimensional sub-algebra 



v 7 = a(x,t) — . (25) 



Assume tt = /( wf+al ' r(i+fl0 * s an exac ^ solution of Eq. (4). Then with the proposed fractional 
method of characteristics, solving the above symmetry equations, we can derive 



v 



(1) = f , xf _ 



r(l+o) °>r(l+a)< 



7,(2) - f( ■£ t a p \ 

,(3) - < 



p etl_xP__ t a \ 

c J vrfi+/3)' r(i+a)/' 

(4)_ f /- a:' 3 -e f Q r -2e\ 
'-/lr(l+/3) e T(l+a) e J 



U 

u 



U 



e •/ vr(i+/3) Zfc r(i+ Q )' r(n-Q)^' 

ffil i z xjg £ T(l+a) 

7/1,0) _ 1 r r(uMirn^U4rf°rn-.')ii 

r ; 73 ° 

V /l+4e r(T+^) 

y « r(l+a)^ T n 

^ l >r(l+/3)r(l+a)+4 £ r(l+a)a;' 3 ' r(l+/3)+4er(l+Q)t a /' 

which are all the classification of solutions of Eq. (4) . 
Take the solution u^ 5 ' as an example, 

„(») =eTO -r™ /( ___ 2e __ i __ ) . (26) 

Assume /( wf+lT) ~^ £ r(i+o:i ' TTT+qy) = c ' wn i c ^ can be set as the initial value of Eq. (4). Now we 

(5) t^e 2 x a s p a /g-, 

can check that u\ = ce r ( 1+/3 ' r < 1+a ) is one of the exact solutions. If we make /( wfr ] 1 rci+ 1 ) = n i 

^e 2 t"e Cg-, 

ce r(i+/3) r(i+oo ; we can derive a new iteration solution u 2 ■ As a result, by similar manipulations, we 
can give u^> . . . u n which are new exact solutions of Eq. (4) . 
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4 Conclusions 

Fractional differential equations have caught considerable attention due to their various applica- 
tions in real physical problems. However, there is no systematic method to derive the exact solution. 
Now, the problem is partly solved in this paper. 

Another problem may arise: can the Lie group method be extended to fractional differential 
equations of fractional order ~ 2? We will discuss such work in future. 

References 

1 K.B. Oldham, J. Spanier, The fractional calculus, Academic Press, New York (1999). 

2 A. A. Kilbas, H.M. Srivastava, J.J. Trujillo, Theory and Applications of Fractional Differential Equations, Elsevier 
Science B.V, Amsterdam, 2006. 

3 V. Kiryakova, Generalized Fractional Calculus and Applications, Longman Scientific & Technical, Harlow, 1994, 
copublished in the United States with John Wiley & Sons, Inc., New York. 

4 V. Lakshmikantham, A.S. Vatsala, Basic theory of fractional differential equations, Nonlinear Anal. 69 (2008) 2677- 
2682. 

5 K.S. Miller, B. Ross, An Introduction to the Fractional Calculus and Differential Equations, John Wiley, New York, 
1993. 

6 I. Podlubny, Fractional Differential Equation, Academic Press, San Diego, 1999. 

7 P.J. Olver, Applications of Lie Groups to Differential Equations, second ed., GTM 107, Springer, Berlin, 1993. 

8 G.W. Bluman, S.C. Anco, Symmetry and integration methods for differential equations, Appl. Math. Sci. Vol. 154, 
Springer, New York, 2002. 

9 E. Buckwar, Y. Luchko, Invariance of a Partial Differential Equation of Fractional Order under the Lie Group of 
Scaling Transformations, J. Math. Anal. Appl. 227 (1998) 81-97. 

10 R.K. Gazizov, A. A. Kasatkin, SY. Lukashchuk, symmetry properties of fractional diffusion equations, Phys. Scr. 
(2009) 014016. 

11 V.D. Djordjevic, T.M. Atanackovic, Similarity solutions to nonlinear heat conduction and Burgers/Korteweg-deVries 
fractional equations, J. Comput. Appl. Math. 222 (2008) 701-714. 

12 H.G. Sun, W. Chen, H. Sheng, Y.Q. Chen, On mean square displacement behaviors of anomalous diffusions with 
variable and random orders, Phys. Lett. A 374 (2010) 906-910. 

13 G. Jumarie, Modified Riemann-Liouville derivative and fractional Taylor series of non-differentiable functions further 
results, Comput. Math. Appl. 51 (2006) 1367-1376. 

14 G. Jumarie, Laplace's transform of fractional order via the Mittag-Leffler function and modified Riemann-Liouville 
derivative, 22 (2009) 1659-1664. 

15 R. Almeida, A.B. Malinowska, D. F. M. Torres, A fractional calculus of variations for multiple integrals with 
application to vibrating string, J. Math. Phys. 51 (2010) 033503. 

16 G.C. Wu, E.W.M. Lee, Fractional Variational Iteration method and Its Appliation, Phys. Lett. A 374 (2010) 2506- 
2509. 

17 G. Jumarie, Lagrange characteristic method for solving a class of nonlinear partial differential equations of fractional 
order, Appl. Math. Lett. 19 (2006) 873-880. 

18 F. Mainardi, Fractional relaxation-oscillation and fractional diffusion-wave phenomena, Chaos. Soliton. Fract. 7 
(1996) 1461-1477. 

19 F. Mainardi, The fundamental solutions for the fractional diffusion-wave equation, Appl. Math. Lett. 9 (1996) 23-28. 

20 O.P. Agrawal, Solution for a fractional diffusion-wave equation defined in a bounded domain, Nonlinear Dynam. 29 
(2002) 145-155. 

21 K. Al-Khaled, S. Momani, An approximate solution for a fractional diffusion- wave equation using the decomposition 
method, Appl. Math. Comput. 165 (2005) 473-483. 

22 N. Ozdemir, D. Karadeniz, Fractional diffusion-wave problem in cylindrical coordinates, Phys. Lett A 372 (2008) 
5968-5972. 

23 S. Das, Analytical Solution Of A Fractional Diffusion Equation By Variational Iteration Method, Comput. Math. 
Appl. 57 (2009) 483-487. 



